Search for: All records

Projecting climate change is a generalization problem: We extrapolate the recent past using physical models across past, present, and future climates. Current climate models require representations of processes that occur at scales smaller than model grid size, which have been the main source of model projection uncertainty. Recent machine learning (ML) algorithms hold promise to improve such process representations but tend to extrapolate poorly to climate regimes that they were not trained on. To get the best of the physical and statistical worlds, we propose a framework, termed “climate-invariant” ML, incorporating knowledge of climate processes into ML algorithms, and show that it can maintain high offline accuracy across a wide range of climate conditions and configurations in three distinct atmospheric models. Our results suggest that explicitly incorporating physical knowledge into data-driven models of Earth system processes can improve their consistency, data efficiency, and generalizability across climate regimes.

 

Abstract Tropical areas with mean upward motion—and as such the zonal-mean intertropical convergence zone (ITCZ)—are projected to contract under global warming. To understand this process, a simple model based on dry static energy and moisture equations is introduced for zonally symmetric overturning driven by sea surface temperature (SST). Processes governing ascent area fraction and zonal mean precipitation are examined for insight into Atmospheric Model Intercomparison Project (AMIP) simulations. Bulk parameters governing radiative feedbacks and moist static energy transport in the simple model are estimated from the AMIP ensemble. Uniform warming in the simple model produces ascent area contraction and precipitation intensification—similar to observations and climate models. Contributing effects include stronger water vapor radiative feedbacks, weaker cloud-radiative feedbacks, stronger convection-circulation feedbacks, and greater poleward moisture export. The simple model identifies parameters consequential for the inter-AMIP-model spread; an ensemble generated by perturbing parameters governing shortwave water vapor feedbacks and gross moist stability changes under warming tracks inter-AMIP-model variations with a correlation coefficient ∼0.46. The simple model also predicts the multimodel mean changes in tropical ascent area and precipitation with reasonable accuracy. Furthermore, the simple model reproduces relationships among ascent area precipitation, ascent strength, and ascent area fraction observed in AMIP models. A substantial portion of the inter-AMIP-model spread is traced to the spread in how moist static energy and vertical velocity profiles change under warming, which in turn impact the gross moist stability in deep convective regions—highlighting the need for observational constraints on these quantities. Significance Statement A large rainband straddles Earth’s tropics. Most, but not all, climate models predict that this rainband will shrink under global warming; a few models predict an expansion of the rainband. To mitigate some of this uncertainty among climate models, we build a simpler model that only contains the essential physics of rainband narrowing. We find several interconnected processes that are important. For climate models, the most important process is the efficiency with which clouds move heat and humidity out of rainy regions. This efficiency varies among climate models and appears to be a primary reason for why climate models do not agree on the rate of rainband narrowing. 

A large spread in model estimates of the equilibrium climate sensitivity (ECS), defined as the global mean near-surface air-temperature increase following a doubling of atmospheric CO₂concentration, leaves us greatly disadvantaged in guiding policy-making for climate change adaptation and mitigation. In this study, we show that the projected ECS in the latest generation of climate models is highly related to seasonal variations of extratropical low-cloud fraction (LCF) in historical simulations. Marked reduction of extratropical LCF from winter to summer is found in models with ECS > 4.75 K, in accordance with the significant reduction of extratropical LCF under a warming climate in these models. In contrast, a pronounced seasonal cycle of extratropical LCF, as supported by satellite observations, is largely absent in models with ECS < 3.3 K. The distinct seasonality in extratropical LCF in climate models is ascribed to their different prevailing cloud regimes governing the extratropical LCF variability.

 

Accurate precipitation monitoring is crucial for understanding climate change and rainfall-driven hazards at a local scale. However, the current suite of monitoring approaches, including weather radar and rain gauges, have different insufficiencies such as low spatial and temporal resolution and difficulty in accurately detecting potentially destructive precipitation events such as hailstorms. In this study, we develop an array-based method to monitor rainfall with seismic nodal stations, offering both high spatial and temporal resolution. We analyze seismic records from 1825 densely spaced, high-frequency seismometers in Oklahoma, and identify signals from nine precipitation events that occurred during the one-month station deployment in 2016. After removing anthropogenic noise and Earth structure response, the obtained precipitation spatial pattern mimics the one from a nearby operational weather radar, while offering higher spatial (~ 300 m) and temporal (< 10 s) resolution. We further show the potential of this approach to monitor hail with joint analysis of seismic intensity and independent precipitation rate measurements, and advocate for coordinated seismological-meteorological field campaign design.

 

Abstract Daily precipitation extremes are projected to intensify with increasing moisture under global warming following the Clausius-Clapeyron (CC) relationship at about $$ 7\% /^\circ {\text{C}} $$ 7 % / ∘ C . However, this increase is not spatially homogeneous. Projections in individual models exhibit regions with substantially larger increases than expected from the CC scaling. Here, we leverage theory and observations of the form of the precipitation probability distribution to substantially improve intermodel agreement in the medium to high precipitation intensity regime, and to interpret projected changes in frequency in the Coupled Model Intercomparison Project Phase 6. Besides particular regions where models consistently display super-CC behavior, we find substantial occurrence of super-CC behavior within a given latitude band when the multi-model average does not require that the models agree point-wise on location within that band. About 13% of the globe and almost 25% of the tropics (30% for tropical land) display increases exceeding 2CC. Over 40% of tropical land points exceed 1.5CC. Risk-ratio analysis shows that even small increases above CC scaling can have disproportionately large effects in the frequency of the most extreme events. Risk due to regional enhancement of precipitation scale increase by dynamical effects must thus be included in vulnerability assessment even if locations are imprecise. 

Conditional instability and the buoyancy of plumes drive moist convection but have a variety of representations in model convective schemes. Vertical thermodynamic structure information from Atmospheric Radiation Measurement (ARM) sites and reanalysis (ERA5), satellite-derived precipitation (TRMM3b42), and diagnostics relevant for plume buoyancy are used to assess climate models. Previous work has shown that CMIP6 models represent moist convective processes more accurately than their CMIP5 counterparts. However, certain biases in convective onset remain pervasive among generations of CMIP modeling efforts. We diagnose these biases in a cohort of nine CMIP6 models with subdaily output, assessing conditional instability in profiles of equivalent potential temperature,θ_e, and saturation equivalent potential temperature,θ_es, in comparison to a plume model with different mixing assumptions. Most models capture qualitative aspects of theθ_esvertical structure, including a substantial decrease with height in the lower free troposphere associated with the entrainment of subsaturated air. We define a “pseudo-entrainment” diagnostic that combines subsaturation and aθ_esmeasure of conditional instability similar to what entrainment would produce under the small-buoyancy approximation. This captures the trade-off between largerθ_eslapse rates (entrainment of dry air) and small subsaturation (permits positive buoyancy despite high entrainment). This pseudo-entrainment diagnostic is also a reasonable indicator of the critical value of integrated buoyancy for precipitation onset. Models with poorθ_e/θ_esstructure (those using variants of the Tiedtke scheme) or low entrainment runs of CAM5, and models with low subsaturation, such as NASA-GISS, lie outside the observational range in this diagnostic.

 

Abstract The performance of GCMs in simulating daily precipitation probability distributions is investigated by comparing 35 CMIP6 models against observational datasets (TRMM-3B42 and GPCP). In these observational datasets, PDFs on wet days follow a power-law range for low and moderate intensities below a characteristic precipitation cutoff scale. Beyond the cutoff scale, the probability drops much faster, hence controlling the size of extremes in a given climate. In the satellite products analyzed, PDFs have no interior peak. Contributions to the first and second moments tend to be single-peaked, implying a single dominant precipitation scale; the relationship to the cutoff scale and log-precipitation coordinate and normalization of frequency density are outlined. Key metrics investigated include the fraction of wet days, PDF power-law exponent, cutoff scale, shape of probability distributions, and number of probability peaks. The simulated power-law exponent and cutoff scale generally fall within observational bounds, although these bounds are large; GPCP systematically displays a smaller exponent and cutoff scale than TRMM-3B42. Most models simulate a more complex PDF shape than these observational datasets, with both PDFs and contributions exhibiting additional peaks in many regions. In most of these instances, one peak can be attributed to large-scale precipitation and the other to convective precipitation. Similar to previous CMIP phases, most models also rain too often and too lightly. These differences in wet-day fraction and PDF shape occur primarily over oceans and may relate to deterministic scales in precipitation parameterizations. It is argued that stochastic parameterizations may contribute to simplifying simulated distributions. 

Abstract A set of diagnostics based on simple, statistical relationships between precipitation and the thermodynamic environment in observations is implemented to assess phase 6 of the Coupled Model Intercomparison Project (CMIP6) model behavior with respect to precipitation. Observational data from the Atmospheric Radiation Measurement (ARM) permanent field observational sites are augmented with satellite observations of precipitation and temperature as an observational baseline. A robust relationship across observational datasets between column water vapor (CWV) and precipitation, in which conditionally averaged precipitation exhibits a sharp pickup at some critical CWV value, provides a useful convective onset diagnostic for climate model comparison. While a few models reproduce an appropriate precipitation pickup, most models begin their pickup at too low CWV and the increase in precipitation with increasing CWV is too weak. Convective transition statistics compiled in column relative humidity (CRH) partially compensate for model temperature biases—although imperfectly since the temperature dependence is more complex than that of column saturation. Significant errors remain in individual models and weak pickups are generally not improved. The conditional-average precipitation as a function of CRH can be decomposed into the product of the probability of raining and mean precipitation during raining times (conditional intensity). The pickup behavior is primarily dependent on the probability of raining near the transition and on the conditional intensity at higher CRH. Most models roughly capture the CRH dependence of these two factors. However, compensating biases often occur: model conditional intensity that is too low at a given CRH is compensated in part by excessive probability of precipitation. 

Abstract In convective quasi-equilibrium theory, tropical tropospheric temperature perturbations are expected to follow vertical profiles constrained by convection, referred to as A-profiles here, often approximated by perturbations of moist adiabats. Differences between an idealized A-profile based on moist-static energy conservation and temperature perturbations derived from entraining and nonentraining parcel computations are modest under convective conditions—deep convection mostly occurs when the lower troposphere is close to saturation, thus minimizing the impact of entrainment on tropospheric temperature. Simple calculations with pseudoadiabatic perturbations about the observed profile thus provide useful baseline A-profiles. The first EOF mode of tropospheric temperature (TEOF1) from the ERA-Interim and AIRS retrievals below the level of neutral buoyancy (LNB) is compared with these A-profiles. The TEOF1 profiles with high LNB, typically above 400 hPa, yield high vertical spatial correlation (∼0.9) with A-profiles, indicating that tropospheric temperature perturbations tend to be consistent with the quasi-equilibrium assumption where the environment is favorable to deep convection. Lower correlation tends to occur in regions with low climatological LNB, less favorable to deep convection. Excluding temperature profiles with low LNB significantly increases the tropical mean vertical spatial correlation. The temperature perturbations near LNB exhibit negative deviations from the A-profiles—the convective cold-top phenomenon—with greater deviation for higher LNB. In regions with lower correlation, the deviation from A-profile shows an S-like shape beneath 600 hPa, usually accompanied by a drier lower troposphere. These findings are robust across a wide range of time scales from daily to monthly, although the vertical spatial correlation and TEOF1 explained variance tend to decrease on short time scales. 

Humid‐heat extremes threaten human health and are increasing in frequency with global warming, so elucidating factors affecting their rate of change is critical. We investigate the role of wet‐bulb temperature (T_W) frequency distribution tail shape on the rate of increase in extremeT_Wthreshold exceedances under 2°C global warming. Results indicate that non‐GaussianT_Wdistribution tails are common worldwide across extensive, spatially coherent regions. More rapid increases in the number of days exceeding the historical 95th percentile are projected in locations with shorter‐than‐Gaussian warm side tails. Asymmetry in the specific humidity distribution, one component ofT_W, is more closely correlated withT_Wtail shape than temperature, suggesting that humidity climatology strongly influences the rate of future changes inT_Wextremes. Short non‐GaussianT_Wwarm tails have notable implications for dangerous humid‐heat in regions where current‐climateT_Wextremes approach human safety limits.